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The growth and assembly of a massive galaxy at z ~ 2 
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ABSTRACT 

We study the stellar mass assembly of the Spiderweb Galaxy (MRC 1138-262), a massive 
z = 2.2 radio galaxy in a protocluster and the probable progenitor of a brightest cluster galaxy. 
Nearby protocluster galaxies are identified and their properties are determined by fitting stellar 
population models to their rest-frame ultraviolet to optical spectral energy distributions. We 
find that within 150 kpc of the radio galaxy the stellar mass is centrally concentrated in the 
radio galaxy, yet most of the dust-uncorrected, instantaneous star formation occurs in the 
surrounding low-mass satellite galaxies. We predict that most of the galaxies within 150kpc 
of the radio galaxy will merge with the central radio galaxy by z — 0, increasing its stellar 
mass by up to a factor of ~ 2. However, it will take several hundred Myr for the first mergers 
to occur, by which time the large star formation rates are likely to have exhausted the gas 
reservoirs in the satellite galaxies. The tidal radii of the satellite galaxies are small, suggesting 
that stars and gas are being stripped and deposited at distances of tens of kpc from the central 
radio galaxy. These stripped stars may become intracluster stars or form an extended stellar 
halo around the radio galaxy, such as those observed around cD galaxies in cluster cores. 

Key words: galaxies: elliptical and lenticular, cD - galaxies: individual: MRC 1138-262 - 
galaxies: evolution - galaxies: high-redshift 



1 INTRODUCTION 

High redshift radio galaxies (Hz RGs) are among the most luminous 
and massive known galaxies dSevmour et alj|2007l and references 
therein). They contain vast reservoirs of atomic and molecular gas, 
and are commonly embedded in giant ionized gas halos. A number 
of well studied HzRGs have been found to be undergoing violent 
activity such as extreme star formation and merging. They contain 
active galactic nuclei (AGN) implying the co-evolution of their su- 
permassive black holes and stellar hosts. Some radio galaxies have 
been found to be located in dense regions of the early Universe that 
may collapse to form clusters of gala xies with masses grea ter than 
10 14 Mq when extrapolated to z = dVenemans et al J 20071) . There- 
fore HzRGs are likely to be the progenitors of low-red s hift bright- 
est cluster galaxies (BCGs; iMilev & De Breuckl 120081 ; IZirm et all 
l2005l:|Pentericci et alll997l) . 

With a mass of ~ 10 12 M W , the radio galaxy MRC 1138- 
262 at z = 2.2 (Fig . 1; IPentericci et alj|l997l ; llVlilev et alj l2006t 
I Overzier et al J 120071) . also known as the ' Spiderweb Galaxy', i s 
one of the most massive galaxies at z ~ 2 ( Sey mour et alj|2007h . 
It lies in a protocluster identified by an overdensity of Lya and 
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Hce emitters, extremely red objects and sub-mm bright galax- 
ies JPentericci et al.ll200(j|: iKurk et al]|2004 iKodama et all 120071 : 
IStevens et alj|2003l : IZirm et alj|2008h . The radio galaxy is 2 mag 
brighter in Ks than any oth er protocluster g alaxy candidate within 
a 1.5 Mpc physical radius dKurk et al] |2004). It is therefore highly 
likely that this radio galaxy is the central galaxy within a massive 
dark matter halo at the centre of the protocluster, and that the nearby 
galaxies may become satellite galaxies, or eventually merge with it 
to form a brightest cluster galaxy. 

Hubble Spac e Telescope (HST) images of this radio source in 
iMilevetal .1 (2006) show that MRC 1 138-262 is surrounded by tens 
of UV-bright clumps and galaxies. These observations agree qual- 
itatively with predictions of hierarchical galaxy formation theory. 
In this work we examine the properties of the nearby galaxies in 
more detail, in order to quantify the future growth and mass assem- 
bly of the radio galaxy through the accretion of these galaxies. We 
combine high resolution rest-frame UV and optical data to identify 
protocluster galaxies that are likely to merge with the radio galaxy. 
By comparing their photometry to stellar population models, we es- 
timate the stellar mass of these galaxies. We then use the observed 
stellar mass distribution as initial conditions to examine the future 
growth and assembly of the radio galaxy. 

In section [2] we describe the observations and data reduction 
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and section[3]describes how the protocluster galaxies are identified. 
Section [4] gives details of the spectral energy distribution (SED) 
fitting, through which we derive the masses of the galaxies. These 
properties are discussed and compared to semi-analytic models in 
section [5] In section [6] we use simple analytic approximations of 
the merging timescale to examine the evolution of MRC 1 138-262. 
Throughout this work we use Hp = 71 kmr 1 , Q.m = 0.27, and 
fl A = 0.73 ( ISpergel et alj|2003h , and at the redshift of MRC 1 138- 
262 (z=2.156) the linear scale is 8.4 kpc arcsec . All magnitudes 
are AB magnitudes. 



2 OBSERVATIONS 

The protocluster surrounding MRC 1138-262 was observed with 
the Advanced camera for Surveys (ACS) instrument on the HST 
though the F475W (g47s) filter for 9 orbits, and 10 orbits though the 
F814W (7gi4) filter. The data w ere reduced using the ACS pipeline 
Blakeslee et alj|2003l) . Details of t he da ta reduction are 



Milev et all d2006t) and lHatch et alj J2008h . The protoclus- 



(APSIS; 
given in 

ter was also observed with NICMOS camera 3 on HST for 3 orbits 
though each of the Fl 10W (/j 10 ) and F160W (H m ) filt ers. Further 
discus sion of the NICMOS observations can be found in lZirm et al.l 
d2008h . 

The HST observations were supplemented with a 3 hour expo- 
sure in U„— band from the LRIS-B instrument on the Keck tele- 
scope (P.I. W. van Breugel), and near-infrared (near-IR) images 
from the Very Large telescope (VLT) ISAAC instrument: Ks— band 
(1.6hour) and a narrow-band image centred on 2.07 \im (4.8hour). 
The U„ —band data has a full width at half maximum (FWHM) of 
1.2 arcsec, and both near-IR images have a FWHM of 0.45 arcsec, 
measured from bright, unsaturated stars that lie close to the radio 
galaxy. Det ails of the reduct ion of the near-IR observations can 
be found in iKurk et al.l j2004l) . Ha emission from the protocluster 
at z= 2.156 lies within the bandwidth of the narrow-band NB207 
filter. Both an Ha emission-line image of the protocluster and an 
emission-line-free Ks continuum image were created from the near- 
IR dataset. 

The images were registered and matched in pixel-scale 
(0.2 arcsec pixel -1 ) and resolution to the Ks image. The 3o im- 
age depths of the images are 27.7 mag in #475, 27.7 mag in 
Ig\4, 26.2 mag in 7no, 26.6 mag in //igo, 24.4 mag in Ks and 
26.8 mag in U„. Magn itudes were measured using SExtractor 
jBertin& Arnoutslll996l) in double-image mode, using a "white- 
light" image as a detection image. The detection image was cre- 
ated by combining the high signal-to-noise g475, I^u, Jno and //jeo 
HST images. Regions with 7 or more connecting pixels which had 
a signal-to-noise ratio greater than 4 times the local background 
noise were identified and isophotal magnitude were measured. 

The point spread function (PSF) of the U„ —band data is more 
than twice as large as the PSF of the other images. To preserve the 
high resolution of the HST and Ks images we scale the U n mag- 
nitudes by a factor to account for the light dispersed beyond the 
galaxy isophotal detection area by the larger PSF. The #475 was 
smoothed with a Gaussian of FWHM 1.2 arcsec to match the res- 
olution of the U„ image. Magnitudes were obtained from both the 
g475 and U„ images using the original detection image and multi- 
plied by the ratio of flux within the isophotal area of the 0.45 arcsec 
resolution £475 image to the 1.2 arcsec resolution 5475 image. Ob- 
jects which are fainter in the 0.45 arcsec resolution image than the 
1.2 arcsec resolution £475 image are blended and no U n magnitude 



was assigned to these objects. All magnitudes were co rrected for 
Galactic extinction of A B = 0. 172 dSchlegel et all 1998b . 



2.1 Emission-line contribution 

Lya emitted from galaxies in the protocluster lies within the band- 
width of the U„ filter, although the filter transmission a t this wave- 
length is only 30% of the maximum filter transmission. IKurk et al.l 
(2004) lists the bright Lya emitters close to the radio galaxy. 
Galaxy 10 (see Fig.[2j has an equivalent width (EW) greater than 
1000A. The U„ magnitude for this galaxy is likely to be dominated 
by the Lya emission. We therefore assume the observed magni- 
tude is an upper limit on the U n continuum emission. The other 
Lya emitters have EWs of 30-100A, therefore a maximum of 5% 
of the emission in the U n — band results from Lya line emission and 
we do not include any correction for this. 

Both the Civ and He II emission lines fall within the £475 
bandwidth, and co ntribute greatly to the diffuse intergalactic light 
(Hatch et al. 2008]), but spectroscopy shows that neither line is 
detected in the majority of the Lya emitting galaxies close to 
MRC 1 138-262 jKuridl2003l) . Only the radio galaxy and galaxy 10 
have large combined Civ and He II EWs of 101A and 87A re- 
spectively and the #475 magnitudes are reduced by 22.5% and 20% 
respectively. 

H/3 and [Olll] emitted from the protocluster galaxies falls 
within the #160 bandwidth. The contribution in H\(,q from these 
lines is estimated using the Ha image and observations of 
MRC 1 138-262 from the integral field unit spectrometer SINFONI, 
which show that Ha+[Nll]> 2.4x[OlIl]+H/3 (N. Nesvadba - pri- 
vate communication). The line emission from the galaxies is likely 
to be fainter than the line emission from the halo, therefore a maxi- 
mum of 5% of the #160 fl ux can be due to line emission. As we do 
not know exactly how much line emission should be removed from 
each galaxy, we do not remove this component. The Ha emission 
was removed from the Ks image using the narrow-band NB207 im- 
age. 



2.2 Removing the nuclear contribution 

The nuclear radio emission is associated wit h galaxy 1 and we refe r 
to this galaxy as the radio galaxy (see Fig. 2; |Pentericci et alj 19971) . 
To remove the point source contribution from the active galactic nu- 
cleus, all galaxies, except for the radio galaxy, are masked out from 
the original high-resolution HST images, and the resulting image is 
fitted with an ellipsoid model to obtain the best-fitting centroid, po- 
sition angle and axial ratio. The azimuthally averaged radial profile 
is fitted with a classical de Vaucouleurs profile plus a point source 
smoothed to the resolution of the image. The centroid of the point 
source is not constrained to lie at the centre of the large-scale ellip- 
soid. The point source contribution to the Ks, H\^q, Jhq and /gi4 
images is 15%, 30%, 23%, 9% respectively and the galaxy profile 
in #475 and Un— bands is consistent with no point source contri- 
bution. Most AGN have nuclear starbursts and therefore a small 
percentage of star formation is likely being removed together with 
the point source. 

The radio galaxy radial profile in Jhq and #160 is well de- 
scribed by a central point source plus a de Vaucouleurs profile with 
a characteristic radius of 0.6 arcsec (5 kpc - comparable to low red- 
shift cD and giant Elliptical galaxies). However, in the rest-frame 
UV light (7gi4 and #475) the radio galaxy is not well described by 
the de Vaucouleurs plus point source profile at large radii because 
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Figure 1. Combined hu-Ji\o-K s image of MRC 1 138-262 (centre) and its 
neighbouring satellites galaxies within 150 kpc. Total extent of the image 
is 300 kpc on each side. Galaxies identified to be within the protocluster 
are marked and labelled in Fig .[2] Note the bridge of red light connecting 
galaxies 17 and 18 possibly arising from tidal stripping. 

of the large contribution of extended diffuse light (see lHatch et al.l 
12008b . 
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Figure 2. Galaxies selected to be at the same redshift as the radio galaxy 
are outlined in black contours. Numbers correspond to galaxy IDs in Table 
[T]and Fig.[3] The Lya halo of the radio galaxy is shown in grey contours. 
The red galaxies with little current star formation (galaxies 2, 4 and 16) 
are located towards the edge of the Lya halo, whilst those with large star 
formation rates are located close to the radio galaxy. The dotted circle en- 
compasses the area within 150 kpc of the centre of the radio nucleus. 



3 SELECTING PROTOCLUSTER GALAXIES 

We wish to study the central region of the protocluster, the radio 
galaxy and the surrounding galaxies that will interact or merge with 
the radio galaxy. To minimize the contamination of more distant 
protocluster galaxies, only galaxies within a radius of ~ 150 kpc 
(18 arcsec) of the radio galaxy were selected. This is the observed 
maximum radial extension of the Lya halo which surrounds the 
radio galaxy and corresponds to the virial radius of a ~ 10 13 M© 
object (see Fig- [2j». We detect 64 separate objects within this area. 

Eight of the 64 objects are confirmed to be protocluster galax- 
ies through the detection of Lya emissi on (galaxies 1, 5, 6, 7, 9, 
10, 11, 15: lKurldl2003| ; iKurk et alj|2004 see Fig.©. Five objects 
are obs erved to em it Ha emission, galaxies 1,10,11,12 are identi- 
fied by IKurk et al.l d2004l) and we identify a further galaxy (ID 8) 
at 11:40:48.37 -26:29:05.2 (J2000) with an EW of 153A. In to- 
tal 9 satellites of the radio galaxy are identified by line emission 
detections. 

Selecting protocluster galaxies through line emission alone bi- 
ases the sample toward star forming galaxies. At the redshift of the 
protocluster the Balmer and 4000A break falls within the J\ \q filter, 
therefore we can ma ke a simple colour cut to select candidate red 
protocluster galaxies. Izirm et aT] d2008| ) show that there is a signif- 
icant overdensity of galaxies in this protocluster field with colours 
1.1 <7no~/fi6o< 2.1. Using the same colour cut we identify 7 red 
galaxies (galaxies 1, 2, 4, 12, 13, 14 and 16) two of which are also 
Ha emitters. 

We include 3 additional galaxies in our catalogue of satel- 
lite galaxies based on their proximity to the radio galaxy and their 
peculiar morphologies. Galaxies 17 and 18 have a bridge of red 
light connecting them which corresponds to a surface brightness 
enhancement in the Lya halo (see Fig[T]and[2} suggesting that they 



are an inte racting system , but the photometric redshift estimation 
code BPZ dBemtezll2000b assigns redshifts of 3.6 and 0.3 respec- 
tively with low probability. Galaxy 19 has a peculiar linear mor- 
phology and lies close in projection to the radio galaxy. A concen- 
tration of chain and tadpole galaxies were found around the radio 
galaxy dMilev et alj r2006) suggesting that they are located in the 
protocluster core. We include galaxies 17, 18 and 19 in the subse- 
quent analysis but note that these galaxies have not been formally 
identified with the protocluster. 

Photometric redshifts of all objects were obt ained using th e 
bayesian photometric redshift estimation code BPZ (Bemtez 2000), 
utilizing the CWWSB template set as described in lBemtezI ( feOOOb 
and 2 additi onal simple stel lar popul ation models with ages of 5 
and 25 Myr dCoe et al.ll2006h from lBruzual & Charlotl d2003b . We 
used a prior based on the galaxy population of the Northern Hubble 
Deep Field and imposed an additional weak prior of a cluster at 
z=2.2 with 1% of the galaxies in the field lying within the cluster. 
This is a weak additional prior since 10% of the galaxies in the 
0.283 arcmin 2 field are identified as protocluster galaxies from line 
emission. 

The Ha emitters labelled 8 and 12, and all protocluster 
galaxies selected by a red colour-cut have photometric redshifts 
within Az=0.3 of the redshift of the radio galaxy. One further 
galaxy (3) has a photometric redshift of 2.2 with Bayesian Odds 
O > 0.9. Including the cluster prior shifted the best fit redshifts of 
galaxies 13, 14 and 16 from z ~ 1.4 to z = 2.2, although z = 2.2 lay 
within the la confidence intervals whether the cluster prior was 
included or not. Increasing the cluster prior such that 50 per cent of 
the galaxies in the field lie within the cluster had no further effect 
and did not increase the number of galaxies with photometric 
redshifts within Az = 0.3 of 2.2. The photometric redshifts of the 
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blue galaxies could not be constrained using BPZ. Galaxies with 
high equivalent width of Lya generally have smaller UV colours 
across the Lyman break than galaxies with no Lya emission. Also, 
although the selection based on U-G is a po werful technique o f 
isolating 'BX' galaxies at z ~ 2.2 ±0.3 (see iReddv et al.l [2005b . 
our U band data does not reach deep enough to constrain the colour 
across the Lyman break accurately. 

Within 150 kpc of the radio galaxy, we find a total of 19 galaxies 
in the protocluster. It should be noted that the n umb er of galax- 
ies see n in the 7gi4 images of iMilev et al.l (2006) and lHatch et al.l 
J2008h is greater than presented in this work. The lower resolution 
of the data used in this work removes many of the smaller clumps 
of bright emission. The large number of the small clumps seen near 
the radio galaxy in the 7gi4 and £475 images suggest they are located 
in the same region as the radio galaxy although no redshift can be 
determined. Thus the 19 detected galaxies should be considered the 
observable subset of a larger population of protocluster galaxies. 



4 STELLAR POPULATION SYNTHESIS MODELING 

A number of works have derived the stellar population properties 
of z ~ 2 galaxies using stellar populatio n synthesis techniques . 
We follow a proced ure s imilar to tha t of Papovich et al.l d200ll) : 
Shaple v et al] d2005h and iKriek et alj i bOoiT and determine stel- 
lar population properties by fitting the broad-band photometry with 
stellar population models. Our aim was to determine the age, intrin- 
sic star formation rate, E(B — V), star formation history and stellar 
mass of each galaxy within 150 kpc of the radio galaxy. 

We compare the spectral energy distribution (SED) of 
the galaxies with stella r population templates created using 
iBruzual & Charlotl d2003h models. We used models with solar 
metallicity, and a Chabrier (2003) initial mass function (IMF) ex- 
tending from 0.1 to 100 Mq. We use the Padova 1994 stellar pop- 
ulation evolution tracks to form single exponentially declining star 
formation histories with characteristic e-folding times (t) of 1, 0.5, 
0.1, 0.05, and 0.01 Gyr, and a model with a constant star-forming 
history. We consider models with 31 ages between 10 6 yrs since 
the onset of star formation and 3.0 Gyr, the age of the Universe at 
z = 2.156, the redshift of MRC 1138-262; and 21 extinction val- 
ues between Ay = and 4 mag, assuming the Calzetti reddening 
law. We compare the photometry to the stellar population synthe- 
sis templates and determine the best-fit parameters by minimum-^ 2 
techniques. Confidence intervals due to the photometric errors are 
derived by constructing 500 synthetic realisations of each galaxy, 
varying the photometry by a random amount drawn from the nor- 
mal photometric error distribution. The stellar population models 
are compared to the synthetic data sets to derive new best-fitting 
models. \a confidence intervals of the best fit parameters are de- 
rived from the minimum and maximum values enclosing 68 per 
cent of the simulations. We also derive best-fit parameters without 
the l/„-band information by comparing the #475 through Ks— band 
photometry to the stellar population models. The results of all the 
alternative fits are within the errors given in Table Q] 

Stellar templates created with single exponentially decreas- 
ing or constant star formation histories provide lower limits to the 
galaxy mass, particularly if the galaxy has a very young age and 
active star formation. Such galaxies have relatively small mass- 
to-light ratios as the observed light is dominated by bright young 
stars, which may hide an older stellar population, that could po- 
tentially be more massive than the brighter young population. To 



estimate the amount of mass that may exist in an older stellar pop- 
ulation we fit the photometry with a two component star forma- 
tion history. One component is fixed to be a 3 Gyr population with 
T = 0.01 Gyr, whilst the other model is fixed to be a 1 Myr old con- 
stant star forming population. The galaxy mass estimated from this 
two-component model is an upper limit on the true stellar mass, lrj 
confidence intervals are drawn from 500 Monte Carlo simulations 
as described above. 

The best-fit parameters depend systematically on the choice 
of metallicity and IMF, for example a Salpeter IMF gives mass-to- 
light ratios that are systematically higher by ~ 1 .5 than the Chabrier 
IMF. Larger metallicities result in redder SEDs, and thus allowing 
the metallicity to decrease, can increase the best-fit extinction, and 
decrease the age of the stellar population. A full discussion of how 
the best-fit ga laxy parameters depen d on the IMF and the metallic- 
ity is given in lPapovich et al.l d200ll) . 

The parameters derived from SED fitting are subject to 
large systematic un certainties and degeneracies and we refer to 
Sha plev et alj [2005) for details. The best-fit age, intrinsic star for- 
mation rate and extinction strongly depend on the assumed star for- 
mation history. The data are unable to constrain the star formation 
histories across our range of parameter space, for all the selected 
galaxies. We therefore infer that the derived age, extinction and in- 
trinsic star formation rate are unreliable, and will not be discussed 
further. In the following discussion we refer only to the derived 
stellar mass (columns 3 and 4 of Table[T](, and the dust-uncorrected 
star formation rate (column 2), determined using the observed £475 
luminosity and 



SFR(M e yr 



-1500 



[ erg s 1 Hz 



(1) 



1.5 x 10 28 

dSalim et al]|2007t) . 

Stellar masses do not suffer the same large systematic 
uncertainties due t o the unconstrained star formation history 
dShaplev et al.ll2005h and th erefore the derived m asses are likely to 
be more robust. We note that Sha plev et alj (2005) use near-infrared 
photometry, which means the derived masses are more reliable. The 
rest-frame V— band is longest wavelength in this work, and there- 
fore we expect larger uncertainties. Although Spitzer mid-infrared 
images are available, the high density of the galaxies in this re- 
gion, and the low resolution of the data, does not allow us to ro- 
bustly measure rest-frame near-infrared magnitudes of the individ- 
ual galaxies. 

Recent work has focused on the contribution to the rest-frame 
near-infrared light of thermally pulsating asymptotic giant branch 
stars, which occurs 0.2-2 Gyrs after the onse t of star formation 
dMaraston et~al]|2006l : IBruzual & Chariol l2003). The reddest band 
we use in this work is equivalent to the rest-frame V, and the contri- 
bution of thermally pulsating asymptotic giant branch stars at this 
wa velength is not important . We have compared our results from 
the lBruzual & Charlotl J2003I) models to the new Chariot & Bruzual 
(2007) models which include the thermally pulsating asymptotic 
giant branch stars. We find no significant change in the results given 
in Table [T] in particular the stellar masses are within the lrj errors 
quoted. 

We note that masses derived from multiband photometry 
are subject to syste matic uncertainties from the assumed IMF 
dConrov et al J 120080 . The uncertainty on the assumed IMF may 
drastically change the derived masses, possibly by more than a fac- 
tor of 10. However, the relative masses of galaxies, which are of 
particular importance to this work, are derived with much greater 
certainty. 
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Figure 3. Broad-band photometry (black filled circles) and lo~ errors are shown for the radio galaxy (ID 1) and 18 satellite galaxies. The point source 
contamination from the AGN has been removed from the radio galaxy. The best-fit stellar population models are shown in grey, the diamond symbols are the 
convolution of the best-fit stellar population models with the filter transmission curves. 



5 RESULTS 

5.1 Properties of the galaxies 

5.1.1 Stellar mass 

The nucleus of the ra dio emission is associated with galaxy 1 
dPentericci et alJ[l997h . The total stellar mass of this galaxy ob- 
tained from a fit to a single stellar population model is 1.1 ±0.2 x 
10 12 Mq which is consistent w ith the upper lim it obtained from the 
rest-frame near-infrared light dSevmour et al.l 120070 . This galaxy 
dominates the central 150 kpc in terms of mass, comprising 85±5 
per cent of the total mass. If one assumes that the mass of the satel- 
lite galaxies is closer to the upper mass limit (listed in column 3 of 
TableQJ the fraction of mass in the radio galaxy reduced to 60±13 
per cent. Rest-frame UV and optical light may be scattered from 
the central active nucleus contaminating the stellar light of galaxy 
1, resulting in an overestimate of the stellar mass. Polarization mea- 



surements are essential to determine the fraction of scattered quasar 
light within each band. The upper mass limit of galaxy 1 is an or- 
der of magnitude larger than previous measurements. We therefore 
always assume in the following discussion, that the mass of this ra- 
dio galaxy is that of the lower mas s limit (1 x lO^M^) which is 
consistent with the measurement of lSevmour et alj d2007l) . 

Galaxy 1 is at least ~10 times more massive than any other 
galaxy within 150 kpc. The sample we have selected is unlikely to 
be a complete inventory of satellite galaxies. From a fit of the ob- 
served Ks mag and the derived upper mass limit, we find that the 
3o" Ks— band detection limit of 24.4 mag corresponds to a mass of 
~6x 10 9 M Q . Our sample is unlikely to be missing any satellite 
galaxies more massive than this mass limit. We identify galaxy 1 as 
the progenitor of the BCG, which already dominates the mass bud- 
get of the protocluster core, and we refer to this galaxy as the radio 
galaxy or the central galaxy. The remaining protocluster galaxies 
within 150 kpc of the radio galaxy are referred to as the satellite 
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ID 


SFR 


Mass 


Mass (upper limit) 
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(M Q yr-') 


(1O 9 M ) 


(1O 9 M ) 


method 


1 
1 


i s f. _i_ n 7 


1 1 nn 4- 9 on 




T \7/v 1-Tr/ RR 
L,ytX,rlU, ±5 ±5 


z 


U.U ± U.4 


,7+13 
12_ 7 


i £,+ 15 


R R 


3 


2.2 i 0.5 


1 1 ± 8 


105_ 54 


Photo-z 


A 

4 


U.U zrz U.4 


t,+24 
31 -16 


34_ 8 


RR 


c 

J 


5.2 i 0.5 
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Table 1. The mass is derived from fitting the photometry to a single expo- 
nentially declining star formation history. The mass upper limit is derived 
from a two-model fit to the photometry, in which one model is maximally 
old, i.e., 3 Gyrs and the other is maximally young (1 Myr). Column 5 lists 
the detection methods by which the galaxy was selected to be in the proto- 
cluster. Lya, Ha are objects which have an excess ofline-emission placing 
them at the same redshift as the radio galaxy, BB indicates galaxies with 
strong Balmer breaks inferred from large observed /no-^160 colours. 

galaxies. The radio and the satellite galaxies are collectively re- 
ferred to as the Spiderweb system. 

All galaxies have sufficient red light to be habouring an older 
stellar population. The maximum masses of the galaxies are gener- 
ally 1-20 times greater than the mass derived from fitting a single 
exponentially declining stellar population. The exception is galaxy 
18 whose SED allows up to 45 times more mass to be hidden be- 
neath the young stellar population that dominates the light in the 
observed UV and optical wavelength range. 

5.1.2 Star formation rates 

The total star formation rate of all the galaxies determined from 
the UV emission is 63±8 Mq yr~' . The star formation is predom- 
inately located within the satellite galaxies rather than in the cen- 
tral galaxy. The radio galaxy contains ~ 25 ± 2% of the star for- 
mation occurring within the whole system. Dust extinction in the 
radio galaxy may skew the fraction of UV light observed in the 
radio galaxy compared to the satellite galaxies. This percentage 
does not take into account the diffus e star format i on tha t lies be- 
tween the radio and satellite galaxies. Hatc h et al] J2008h find that 
the amount of intracluster star formation is similar to the amount of 
star formation that occurs within the galaxies in a similar sized area 
(15 x 18arcsec 2 ). Thus most of the star formation occurs beyond 
the region we have defined as the central galaxy. Only 3 satellite 
galaxies have no detected Far-UV flux, the majority are in a phase 
of rapid star formation. From the current data we cannot determine 
whether galaxies 2, 4, and 13 are passively evolving red galaxies, 
or dusty star-forming galaxies. 

Large amounts of dust-obscured star formation is confirmed 
by the bright sub-millimeter (sub-mm) emission observed in 



this galaxy JStevens et al.l l2003h . with a 850[J.m flux density of 
6.7±2.4 mjy. Thus the total star formation rate of the region close 
to the radio galaxy is potentially much higher than given by the 
dust-uncorrected UV luminosity. The galaxies have red UV colours 
(see Fig. 5 in section I5.2.2t implying extinctions in the approxi- 
mate range 0.1 < E(B — V)< 0.5 depending on the age, star for- 
mation history and metallicity of the stellar population. Correcting 
the UV luminosity for this amount of dust extinction would bring 
the star formation rates derived from the observed UV and sub-mm 
emission into agreement. The 850p.m emission from this location 
is spa tially elongated along a position angle of 72° dStevens et al.l 
2003). The satellite galaxies with extreme star formation and large 
amounts of extinction (galaxies 5, 6, 7, 9, 10, 11 and 15) are elon- 
gated approximately East-West along the sky stretching 15 arcsec 
from end-to-end. Therefore the extended sub-mm emission results 
from source confusion and is the sum of the combined star for- 
mation occurring in the radio and satellite galaxies. High spatial 
resolution sub-mm observations have deconvolved the extended 
sub-mm emission from the high redshift galaxy 4c40.7 into two 
separate compone nts indicating merger induced star formation 
dlvison et al]|2008t) . These optical and near-infrared observations 
show the star formation occurs in multiple compact objects and 
high resolution sub-mm may be able separate the sub-mm emis- 
sion into the various galaxy components. 

5.1.3 Spatial distribution of the satellite galaxies 

Almost all of the satellite galaxies lie within the large Lya halo 
that surrounds the radio galaxy, even those galaxies with little or 
no star formation. Since our observations show a single snapshot 
in time, most of the satellite galaxies' orbits must lie within the 
Lya halo. This gives an approximate orbit radius of 100 kpc. The 
asymmetric Lya halo extension to the East also contains satellite 
galaxies. This strong spatial correlation implies a connection be- 
tween the Lya halo and the satellite galaxies. The Lya halo may 
originate from gas stripped from the satellites, or the galaxies may 
excite the gas, possibly by shocks from their movement through an 
existing gas reservoir, or UV emission from hot young stars. The 
total Lya luminosity of the halo is 2.6x 10 45 erg s~' and would re- 
quire a star formation rate of ~2500M P yr~' to ionize it. The UV 
star formation rates are more than an order of magnitude too small 
to ionize the halo, although the bright sub-mm emission suggests 
a large amount of obscured star formation. It is possible that the 
Lya halo is ionized by the UV emission from young hot stars pro- 
duced during this phase of extreme star formation, or is a mixture 
of both stellar and AGN photoionization I Villar-Marti n et al.l 1 9971 : 
IVillar-Martm et al.ll2007ri . 

5.2 Comparison of the Spiderweb system to simulations 

We com pare our results to two d ifferent semi analytic models 
(SAMs): IPeLucia & Blaizotll2007l (DLB07) and lFont et ai]|2008l 
(F08) which are based on the dark mater only Millennium Simula- 
tionU. We construct catalogues from the SAMs that most resemble 
our own selection criteria. We select all galaxies within lOMpc of 
the most massive (stellar mass) 10 galaxies at z = 2.2, then make a 
further cut selecting only those galaxies that lie within a projected 
distance of 150 kpc from the massive galaxy. We take an average of 
3 orthogonal projections of each of the 10 galaxies. 

1 http://www.mpa-garching.mpg.de/millennium/ 
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Figure 4. Comparison between the derived stellar masses of the observed 
galaxies within 150 kpc of the Spiderweb galaxy (black histograms) and 
the masses of galaxies within the same projected distance from the 10 most 
massive galaxies within 2 SAMs (F08 - blue, and DBL07 - red). Poisson 
errors for the SAMs are given by ~(number in each bin) 1 / 2 x 0.18. Only 
those simulated galaxies which have Mh less that the 3o~ detection limit 
(26.6 mag) are plotted. Derived upper mass limits are plotted as black solid 
circles, and lower mass limits for the satellites are plotted as open circles. 

We compile catalogues of the SDSS g', SDSS (', H and K mag- 
nitudes for all simulated galaxies. Although the filter curves differ 
slightly to the observed #475, /814, #160 and K s filters, this will not 
affect our comparisons as we only consider large (1 mag) bins when 
comparing galaxy magnitudes. We convert the SDSSg' luminosi- 
ties to star formation rates using equation Q] for comparison to the 
observed UV star formation rates. The SAMs predict a large num- 
ber of low luminosity satellite galaxies that may bias any compar- 
ison to the data, therefore we make a luminosity cut equivalent to 
the 3(7 detection limit in the Hi o- 

The average number of galaxies within each projected box is 
27 for the F08 models and 30 for DBL07 compared to the 19 galax- 
ies observed within 150 kpc of MRC 1138-262. The SAMs predict 
too many satellite galaxies above our detection threshold compared 
to our observations. This discrepancy is further compounded by the 
fact that the dark matter halo of MRC 1 138-262 is probably more 
massive than its simulated counterparts in the SAMs. The missing 
simulated satellites are faint and have low masses (see FigQ. Such 
faint, low-mass galaxies may not be resolved in the smoothed im- 
ages, especially if they are not compact. A larger num ber of galax- 
ies are seen in the higher resolution HST ' l%n images of Mile v et aU 
fcOOrj) of the sam e region, and the intracluster light discussed in 
Hatch et al. ( 2008) may be partly due to unresolved, low mass satel- 
lites. It is therefore likely that the observations presented here miss 
some low-mass satellite galaxies because they cannot be resolved, 
and thus the discrepancy is most likely due to missing satellite 
galaxies in the observations, rather than the simulations predicting 
too many satellite galaxies. 



5. 2. 1 Stellar Masses 

The stellar mass of the most massive semi analytic galaxy at = 2.2 
is 2.2 x 10 11 Mq in the F08 models and 2.8 x 10 11 M Q in the 
DLB07 models. These are ~4 times less massive than the radio 
galaxy (galaxy 1). Our derived mass of 10 12 M,r ; for the radio 
galaxy places it as one of the most massive galaxies detected at 
z > 2. However, we note that this mass is likely to be an upper 



limit, as we have not taken into account light scattered from the 
central AGN. Mass derived from fitting stellar templates to broad- 
band photometry are systematically overestimated by a factor of 
1.4 — 2 compare d to masses derived from fitting stell ar templates to 
spectroscopy JShaplev et all 2005; Kriek et al. 2008). We must also 
consider that stellar population synthesis modeling contains large 
systematic uncertainties and the mass of these galaxies, especially 
galaxy 1, may be lower than the presented in Table [T] Most high 
redshift radio galaxies at z ~ 2 have similar masses to the most 
massive galaxies in the SAMs ^Seymour et al.ll2o"o"7r) . the Spider- 
web Galaxy has an exceptionally high mass. 

In general the central galaxy in each model contains between 
5 and 35 per cent of the stellar mass within 150 kpc of the massive 
galaxy. If we assume the upper mass limits for the satellite galaxies, 
the central galaxy in the Spiderweb system contains 60± 1 3 per cent 
of the stellar mass. Although the stellar mass distribution of the 
SAMs is generally not as centrally concentrated as the Spiderweb 
system, one of the massive galaxies in the F08 SAM contains 35 
per cent of the selected stellar mass within the central galaxy, and is 
therefore compatible with the lower limit allowed by observations. 
Thus the Spiderweb system may not typical, but is compatible with 
of the stellar mass d istribution of ma ssive galaxies in these SAMs. 

The SAMs of ISaro et al.l 120081 based on a dark matter plus 
non-radiative gas particles produce BCGs that are 25 per cent larger 
than those from dark-matter only simulations, and their dominant 
BCGs contains twice as much stellar mass as any satellite galaxy. 
These simulations appear to match the derived mass distribution of 
the Spiderweb system better than the SAMs based on dark-matter 
only. 

Fig.[4]compares the derived upper and lower limits of the stel- 
lar masses of the galaxies within 150 kpc of the radio galaxy to 
those galaxies within 150kpc of the 10 most massive galaxies in 
the SAMs at z = 2.2. Filled circles plot the upper limit masses, 
whilst open circles plot the lower limit masses. We only plot simu- 
lated galaxies with rriH < 26.6, the 3o" detection limit in #i60- 

The stellar masses of the simulated central galaxies are at least 
a factor of 4 less that the observed radio galaxy, therefore we are 
comparing the observations to the masses of galaxies in lower mass 
halos that of the Spiderweb Galaxy. We thus expect the true distri- 
bution of galaxy masses to be shifted to higher masses compared to 
these simulations. Above 10 10 Mq, the SAMs agree with the de- 
rived lower mass limits, whilst the upper mass limits lie to the right 
of the simulated galaxy mass distribution. This is entirely consis- 
tent with our expectations and imply that the derived upper and 
lower mass limits are likely to bracket the true stellar mass of the 
galaxies. Below 10 10 Mq the simulations predict a larger number 
of satellites than are observed. Comparing our upper mass limits 
to the observed Ks magnitude, we find we are only complete to 
masses above ~ 10 9,8 Mq, thus we are likely to be missing many 
of the low-mass satellites with red colours. 

5.2.2 Galaxy colours 

The right-most panel of Fig.[5] shows a comparison of the range of 
UV continuum slope, /3, for the observed and simulated galaxies 
within 150 kpc of the radio galaxy. The simulated galaxies have 
been subjected to the same select criteria as the observations - the 
galaxies must be observed in both the i and g' bands - thus the 
galaxies must be actively forming stars. The red UV colours sug- 
gest significant amounts of dust between us and the galaxies. The 
amount of extinction, E(B — V), that can cause these red colours can 
range from 0. 1 to 0.5 depending on the age and metallicity of the 
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Figure 5. K and G band luminosity functions and UV slope for the observed galaxies (black) and 2 versions of SAMs (F08 - blue, and DBL07 - red). 
Poisson errors for the SAMs are given by ~(number in each bin) 1 / 2 x 0.18. Both luminosity functions stop at the 3o~ detection limit in the considered band. 
The observed K luminosity function matches reasonably well with the models, suggesting that the mass of the simulated galaxies matches the observations. 
However there is a mismatch in the observed #475 and G luminosity functions implying the star formation rates in the simulated galaxies are too low. Both 
models produce galaxy UV slopes that are too blue, although the DBL07 models are a better match to the observations. This implies the simulated galaxies 
which have active star formation have either too little dust or the stellar ages are too low. 



galaxy's stellar population. DLB07 models provide the closest fit to 
the data, however, the simulated galaxies with active star formation 
have too little dust or the stellar ages are too low. 

5.2.3 Star formation rates 

The SAMs predict total star formation rates of 30-70 M yr~ 1 
within 150kpc radius, similar to the observed UV star formation 
rate of 63±8M^,yr (dust-uncorrected) in the Spiderweb sys- 
tem. The galaxy colours show that the dust content of the simulated 
galaxies is less than that of the observed galaxies, and therefore a 
comparison of the intrinsic star formation rates may result in a dis- 
crepancy. 

We make a more detailed comparison by comparing the ob- 
served K and g magnitudes with those from the SAMs in Fig.[5] The 
left panel shows there is general agreement between the observed 
Ks luminosity function and the SAMs. The Ks band is the reddest 
observational band and covers rest-frame V— band light. From the 
selection of observational data used in this work, it is the best tracer 
of galaxy mass, however there can be significant contamination of 
light from young stars. 

The g— band covers rest-frame 1500A light, which can be con- 
verted to star formation rate through equation[T] This comparison is 
shown in the middle panel of Fig.|5] We note that the observed g 4 75 
luminosity function does not have the usual shape of a UV lumi- 
nosity function. Therefore we do not expect the simu lations, which 
repro duce the general galaxy population at z ~ 2 dGuo & White! 
2008), to match the observed rest-frame UV luminosity function. 
Indeed the general shape of the simulated UV luminosity func- 
tion agrees with studies of the z ~ 2.2 luminosity function (e.g., 
ISawicki & Thompson! l200rj) . However this means the simulated 
satellite galaxies typically have lower g magnitudes, and therefore 
lower star formation rates, than the observed satellite galaxies in 
the Spiderweb system. Since the number of simulated galaxies at a 



given mass approximately matches the observations, this discrep- 
ancy is caused by the low star formation rates of the simulated 
galaxies, rather than a lack of satellite galaxies. The SAMs pre- 
dict centrally star formation, whilst the observations show that most 
of the star formation occurs beyond the central radio galaxy. The 
SAMs reproduce the general galaxy population at z ~ 2, but are 
unable to match some of the details of this special system. 

5.2.4 Smooth-particle hydrodynamical simulations 

ISaro et af. (2009) simulate two protoclusters using smooth parti- 
cle hydrodynamical simulations and have compared their simu- 
lations to observations of the Spiderweb Galaxy, and the galax- 
ies within 75kpc of this radio galaxy. The BCGs from these 
smooth-particle hydrodynamical simulations have stellar masses of 
5 x 10 12 /j _1 M and 8 x 10 11 lT x M Q at z = 2. However, these sim- 
ulations do not include AGN feedback, which will act to reduce 
these masses, thus these masses should be considered upper limits. 
The most massive satellite galaxies within the simulations are a fac- 
tor of 4-10 less massive that the already dominant BCG. Therefore 
the stellar mass distributions of the smooth-particle hydrodynami- 
cal simulations are different to that of the SAMs, and they provide 
a better match to the observed stellar mass distribution of the Spi- 
derweb Galaxy. 

The smooth particle hydrodynamical simulations have total 
star formations rates of 600M,^,yr~ 1 and 1750 Mq yr~' within a 
smaller region (within a box of 150h~'kpc on a side) than that 
considered in the data presented here. These high rates are much 
greater than observed in the UV, although the bright sub-mm emis- 
sion from this source (Stev ens et al . 2003) indicates a large amount 
of obscured star f ormation of the or der of a thousand solar masses 
per year. Thus the lSaroet all [2009) results are consistent with the 
data. 

A comparison to the luminosities of the simulated galaxies of 
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ISaro et al.l ( 120091) is not possible, but we note that they match the 
^814 galaxy luminosity function of the galaxies within 75kpc of 
MRC 1138-262 once they include a prescription for dust attenua- 
tion. 



6 DISCUSSION 

6.1 The state of play at z = 2.2 

Within a 150 kpc radius, the high redshift radio galaxy MRC 1 138- 
262 at z ~ 2.2 is surrounded by 18 satellite galaxies with masses 
between 10 s — 10 11 M s . Most of the galaxies within 150 kpc of 
the radio galaxy are UV-bright line-emitters and have large star 
formation rates. These star forming galaxies have similar upper 
limit masses to the galaxies selected through their strong Balmer 
breaks. Thus the galaxies with sufficient gas reservoirs to fuel 
large star formation rates do not necessarily have low stellar 
mas ses. It is known that star fo rmation can be triggered by merging 
(e.g. lMihos & Hernquistlll994l) and it is likely that this is occurring 
on a large scale in this system, with most (> 75 ± 2 per cent) of 
the stars forming in the satellite galaxies, rather than in the central 
galaxy. 

With a mass of 10 12 Mq the central galaxy is already of a 
similar mass to cD galaxies at the centre of nearby galaxy clus- 
ters. The fraction of mass in the Spiderweb system that is con- 
tained within the satellite galaxies, depends on whether we assume 
the stellar masses are better estimated by the lower or upper mass 
limits, given in columns 3 and 4 of Table [TJ respectively. Assum- 
ing upper (lower) mass limits, the combined mass within all the 
satellite galaxies comprises 40 ± 13 (15 ±5) per cent of the total 
stellar mass. The radio galaxy still has room for further growth, 
and may undergo a number of major-mergers before it reaches the 
limit of the most massive, local cD galaxies. However, galaxy 1 
has a Ks— band magnitude that is 2 magnitudes brighter than any 
other protocluster galaxy candidate known within a 1.5 Mpc radius 
dKurk et alj|2004h . and therefore it is unlikely that the protocluster 
cont ains more galaxies of a sim ilar mass at the observed redshift. 

|Pe Lucia & Blaizotl d2007h find that ~ 70 per cent of BCGs at 
z — 2 do not become BCGs at z = 0, but rather end up as the central 
galaxies of massive halos which are less massive than cluster halos. 
The most massive halos at z = 2 are generally not the progenitors 
of z = clusters. Therefore, whilst MRC 1138-262 appears to the 
the brightest protocluster galaxy at z = 2, it is uncertain whether the 
Spiderweb system will evolve into a local BCG. 



6.2 The evolution of the Spiderweb System 

6.2.1 Merging timescales 

A complete N-body numerical investigation of this system is re- 
quired to study the evolution of the Spiderweb system, which com- 
prises of the radio galaxy and all galaxies within a 150 kpc pro- 
jected radius. The data given in Table [TJ and Fig.|2] can be taken 
as initial conditions for such a simulation. However, if we assume 
that the satellite galaxies all lie on circular orbits around the radio 
galaxy, with radii given by their projected radii, we can use analytic 
approximations to gain an idea of the evolution of this system. 

When an orbital point mass (such as the satellite galaxy) trav- 
els through an uniform background mass distribution (such as the 
dark matter halo of the central galaxy), i t experiences a stea dy de- 
celeration known as dynamical friction lChandrasekel ll943l) . The 



merging timescale (T mer ge) of these galaxies can be taken as the 
dynamical friction timescale given by 

x =lA1 Z J^(^±\ (2 ) 
merge InA V M sateUite J 

where the Coulomb logarithm InA = ln(l + Af centra i/M sate nj te ), 
^central is the mass of the central galaxy, M sate nite is the mass of 
the satellite galaxies and the orbital timescale is given by idyn — 
{-rS — )'/ 2 , where R is the distance between the satellite and cen- 
tral galaxy (Bin nev & Tremaind ll987). M centra i/M sate nj te is the ra- 
tio of dark matter halo mass. Matching the abundances of galaxies 
above a given stellar mass to halos above a given dark-matter mass 
is a powerful t echnique employed to relat e the galactic stellar mass 
to a halo mass dConrov & Wechsler 2008), however, this technique 
is not reliable at high redshift (z > 1). Instead we assume that the 
amount of dark matter within each galaxy is directly proportional 
to the galaxy's stellar mass, therefore the ratio M centra i/M sate uj te 
is the ratio of stellar mass. I f the abundance-matching estimates 
of IConrov & Wechslej ( 120081) are correct, the merging timescales 
will be reduced, especially for low mass satellites. M sate ui te is the 
mass of the satellite galaxies and M centra i = 1 x 10 12 Mq, the lower 
limit of the stellar mass of galaxy 1. The total (dark matter plus 
stellar) mass of the central galaxy is estimated from the available 
dynamical information. Spectroscopic observations of galaxies 5, 
7, 9, 11, 15 show they have line-of-sight velocities of 1030 kms , 
210kms~' , 570k ms~ 1 .llSOkms" 1 and 1240 kms -1 relative to 
the central galaxy dMilev et"al1l2006l) which are consistent with a 
cent ral mass of Mj nta ] = 10 13 Mq . 

iBovlan-Kolchin et alj ( 120081) used numerical simulations to 
find a more realistic merging timescale for minor mergers. The 
merging timescale is increased relative to equation|2]due to the sig- 
nificant mass loss of the satellites as they orbit the central galaxy. 
The merging timescale is given by 

W = 0.216^ f^iy' 3 exp[1.9(l- e 2 ) 1 / 2 ], (3) 
InA \ Az sa tellite/ 

where e is the orbital eccentricity. We assume all satellites have cir- 
cular orbits (e = 0). Lower angular momentum orbits have shorter 
merging times, if the satellites have elliptical orbits (0 < e < 1), 
their merging times could decrease by a factor of ~ 7. 

The merging times for all satellite galaxies within 150 kpc of 
the central radio galaxy are shown in Fig.[6] The filled circles mark 
the merging timescales derived using equation[3]and the upper limit 
of the galaxy mass, whilst the open circles use the lower limit of 
the galaxy mass. The central galaxy is always assumed to have a 
mass of 1.1 x 10 12 M e , which is consistent with other measure- 
ments from the rest-frame near-infrared. The open circles indicate 
the longest possible merging times, as we have assumed both the 
lowest masses and orbital eccentricities of the satellite galaxies. 
The vertical lines show the time from z = 2.2 until z= 1 and z= 0. 

The more massive galaxies generally merge earliest, irrespec- 
tive of their initial radial separation from the central radio galaxy. 
Most of the satellite galaxies will merge by z = if they have 
masses close to their upper mass limits. In doing so the radio galaxy 
will increase its mass by a factor of 1.5^5 1. Alternatively, the mass 
of the satellite galaxies may be closer to the lower mass limits. In 
this case, only a third of the satellites will merge before z = 0. Thus 
the amount of growth is dramatically reduced, with only a ~13 per 
cent increase in the stellar mass through merging. However, this is 
not the full story, as tidal stripping is an important process that also 
depends on the mass of the satellite galaxies. The minimum merg- 
ing timescale is ~ 300 Myr which allows the satellites to evolve 
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Figure 6. The merging times for all satellite galaxies within 150kpc of the 
central radio galaxy. The dotted and dashed lines indicate the time from z = 
2.2 until z = 1 and z = respectively. Merging timescales are derived from 
equation[3] Filled circles indicate merging times using the upper mass limits 
from column 4 of Table [T] whilst the open circles are merging timescales 
derived assuming the satellites have the lower masses from column 3. Red 
symbols are galaxies who are likely to be experiencing rapid tidal stripping 
as their tidal radii are less than their half-light diameters. Most satellites 
will either merge with the central radio galaxy by z = 0, or are experiencing 
rapid tidal stripping, and are therefore unlikely to survive until z — as 
separate entities. 

before they merge. Thus many of them are likely to be older, red- 
der and have lower star formation rates by the time they merge. The 
gas reservoirs in the satellites may no longer exist if the satellites 
continue to form stars rapidly. 

6.2.2 Tidal stripping 

The typical orbital time is lOOMyrs, therefore the satellite galax- 
ies will orbit the central galaxy between several and hundreds of 
times before they merge. As the satellites orbit the central galaxy, 
tidal interactions can strip large quantities of gas and stars from the 
satellite galaxies into tidal tails, which are only partly re-acquired 
by the galaxies ddi Matteo et al .1 12007). After several orbits kine- 
matically distinct streams will appear. One possible such stream 
lies between galaxies 17 and 18 (see Figs .Q] and O. 

The tidal radius of the satellites can be approximated by 

R tidcU ~ ( ^satellite \ 3 R (4) 
V ^central ) 

where R is the di stance between the satellit e galaxy and the cen- 
tral radio galaxy l lBinnev & Tremaine1ll987h . This radius approxi- 
mately marks the limit within which stars are bound to the satellite, 
and beyond which stars can be tidally stripped. Fig.|7] shows that 
the tidal radius for many of the satellite galaxies is less than 10 kpc, 
hence the outer stars of the satellites are likely being stripped away. 
Further stripping may occur from interactions between the satel- 
lite galaxies, including those galaxies below the observational de- 
tection limit. The mass ratios of the satellite galaxies are closer to 
unity, which increases the tidal radii. However, the radial separation 
between two satellites may become small during some part of their 
orbits, and during this period the forces between the satellite galax- 
ies may dominate, increasing the amount of tidal stripping that oc- 



Figure 7. The tidal radius of the satellite galaxies derived from the upper 
mass limit of the satellites shown as filled circles, whilst radii derived from 
lower mass limits are open circles. Red points are galaxies whose tidal radii 
are smaller than their half-light diameters as measured from the //i60 image. 
All of these satellites are likely to be experiencing some tidal stripping, but 
those coloured red will be stripped at a faster rate. 



curs. Such interactions can only be modeled by a N-body numerical 
simulation of the system. 

The red symbols in Figs.[6]and[7]mark the galaxies which have 
smaller tidal radii than half-light diameters, measured from the 
#160 image. These galaxies are probably experiencing tidal strip- 
ping of their central stars, and it is unlikely that these galaxies will 
survive as distinct entities for a large number of orbits. It is likely 
that most of their stars will be stripped before they merge. There 
is a strong dependance on mass, so the satellite galaxies would ex- 
perience more stripping if their mass was closer to the lower mass 
limit rather than the upper mass limit. 

Fig.[7] shows that the stars and gas can be stripped from the 
satellite galaxies at large distances, more than 50 kpc from the cen- 
tral galaxy. Therefore these stripped stars may form an extended 
stellar halo around the radio galaxy, so the radio galaxy can grow 
in size, or form something akin to a cD halo. It has been pos- 
tulated that the extended stellar halo s of cD galaxies may form 
as stars are stripped during m ergers dGallagher & Ostriken[l972l 
iMalum ufh & Ri chstonel 1 19841) and these observations may have 
captured the Spiderweb Galaxy in the early stages of such a pro- 
cess. A numerical N-body simulation of this system will determine 
whether such conditions can result in an extended stellar halo such 
as a cD halo. 



6.2.3 Assembly scenarios of the Spiderweb system 

The evolution of the system strongly depends on the mass ratio 
of the central and satellite galaxies. Since we have assumed that 
the central galaxy has a mass of 10 12 M Q , consistent with previous 
measurements, the evolution of this system depends on the mass of 
the satellite galaxies. If the galaxies have true masses similar to the 
derived upper mass limits, most of the satellites will merge by z = 0, 
with the most massive systems merging within a few hundred Myrs. 
Tidal stripping is most important for the low mass satellites that are 
within 30 kpc of the central galaxy. The stars accumulated from the 
mergers and stripping are likely to be deposited close to the central 



© RAS, MNRAS 000, 000-000 



The growth and assembly of a massive galaxy 1 1 



galaxy, and the central galaxy can almost double its z = 2.2 stellar 
mass through these minor mergers. 

Alternatively, if the masses of the satellite galaxies are similar 
to the lower limits we have presented, only the third most massive 
satellites will merge with the central galaxy by z = 0, and most 
merging occurs at later times, between 1 < z < 0. The mass of 
the central galaxy is not greatly increased by these mergers, al- 
though this is due to the low masses of the satellites, rather than 
the low number of mergers. The other satellites are unlikely to sur- 
vive as distinct entities until z = 0. Fig.[6]shows that those satellites 
that will not merge by z = are experiencing strong tidal strip- 
ping of their inner stars. These stars will be stripped from the satel- 
lites and become bound to the central galaxy. In such a manner 
the central galaxy may accumulate all the stellar mass that is cur- 
rently observed in the satellites. The satellites experiencing strong 
tidal stripping generally have the lowest mass, and longest merging 
timescale, and can be stripped at large separations from the central 
galaxy, up to 60 kpc. 

In both scenarios the central galaxy is likely to accumulate 
most of the stellar mass of the satellite galaxies by z = 0, however, 
it does so in different manners. If the satellites' masses are large, 
most of the stellar mass from these satellites is deposited relatively 
quickly and at small radial distances from the core of the central 
radio galaxy. On the other hand, if the satellites have low masses, 
most mass is deposited through relatively slow mergers. A fraction 
of the stellar mass will be deposited at large radii from the core of 
the central galaxy by tidal stripping, forming an extended stellar 
halo. The total growth of the central galaxy due to the consump- 
tion of the satellite galaxies is between -~10 and 120 per cent. This 
large range is due to the large uncertainty of the stellar mass derived 
through stellar template synthesis modeling. The future evolution 
of this system is likely to be somewhere in between the two sce- 
narios present above. In addition, this amount of growth should be 
considered a lower limit, as massive nearby galaxies (> 10 11 M Q ), 
that lie beyond 150 kpc, may also merge with the central galaxy. 

The crucial parameter that drives the evolution of such a sys- 
tem is the mass ratio of the central and satellite galaxies. If a galaxy 
is able to assemble quickly such that it is many times more massive 
than the general galaxy population at high redshift (z > 2), a larger 
fraction of the stellar mass from its future mergers and close en- 
counters will end up in an extended stellar halo. Thus whether or 
not a brightest cluster galaxy becomes a cD galaxy, may depend 
on its early-time formation history. The Spiderweb galaxy, with its 
large mass at z = 2.2, is in a good position to build a substantial 
extended stellar halo by z = 0. 



7 CONCLUDING REMARKS 

We have used multi-wavelength high-resolution images to study 
the assembling radio galaxy, MRC1 138-262, at z ~ 2.2. We iden- 
tified a galaxy population that lies at the same redshift and within 
a projected distance of 150 kpc from the radio galaxy. Fitting stel- 
lar population models to the photometry of these galaxies we have 
examined the current distribution of stellar mass and ongoing star 
formation, and used analytic approximations to predict the future 
evolution of this galaxy. 

1. We identified 19 protocluster galaxies, comprising of a cen- 
tral massive galaxy (10 12 Mq) that is surrounded by 18 smaller 
satellite galaxies with masses between 10 8 - 10 11 M Q . The cen- 
tral galaxy is already of a similar mass to local cD galaxies, but 
there is still room for growth as cD galaxies can have masses of 



10 Mpj. The satellite galaxies contribute between 10 and 50 per 
cent of the stellar mass of the whole system, and they dominate the 
star formation rate, containing 75 ± 2 per cent of the instantaneous 
star formation (dust-uncorrected). 

2. The majority of the satellite galaxies have large star forma- 
tion rates and are inferred to be gas-rich. The intrinsic star forma- 
tion rates are uncertain because the red UV continuum slopes indi- 
cate dust obscuration. The high star formation rates of the satellite 
galaxies coupled with the presence of dust explai ns the spatially ex - 
tended 850^.m emission from MRC 1 138-262 dstevens et al.l2003l) . 
The extended sub-mm emission most likely results from the sum 
of the compact star formation occurring in the radio and satellite 
galaxies rather than a system-wide starburst. These results predict 
that a large fraction of the sub-mm flux should come from the west 
of the radio galaxy and high resolution sub-mm observations may 
be able to separate the emission into the various galaxy compo- 
nents. 

3. We have compared the properties of the galaxies within 
150 kpc of MRC 1138-262 to predictions of semi-analytic models. 
The models are able to reproduce the observed K— band galaxy lu- 
minosity function, and the derived stellar masses of the satellite 
galaxies. The models predict galaxy colours which are bluer than 
the observations, and yet they do not match the observed high star 
formation rates. There may be a process, not accounted for in the 
models, which increases star formation when galaxies are situated 
in a dense merging environment, such as the Spiderweb system. 

4. Most of the satellite galaxies within 150 kpc will accrete 
onto the central radio galaxy before z = 0, increasing its mass by 
between ~10 and 120 per cent. The merging timescales are long, 
and the vast majority of the galaxies will not merge for at least 
300Myrs, by which time they will have completed many orbits. 
Gas is rapidly converted to stars as inferred from the large star for- 
mation rates. Therefore the final merging galaxies are likely to con- 
tain very little gas, and the mergers are likely to be dissipationless. 

5. The tidal radii of the satellite galaxies are small, generally 
less than 10 kpc, and in some cases are less than the half-light di- 
ameter of the galaxies. Therefore stars and gas are probably be- 
ing stripped from the satellite galaxies. This stripping can occur at 
large distances from the central galaxy (more than 60 kpc) and the 
stripped stars may form an extended stellar halo such as those that 
encompass local cD galaxies. 
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